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Abstract 

This work finds the non isotropic noncentral elliptical shape distributions via SVD 
decomposition in the context of zonal polynomials, avoiding the invariant polyno- 
mials and the open problems for their computation. The new shape distributions 
are easily computable and then the inference procedure is based on exact densities 
instead of the published approximations and asymptotic densities of isotropic mod- 
els. An application of the technique is illustrated with a classical landmark data 
in Biology, for this, three models are proposed, the usual Gaussian and two non 
Gaussian; the best one is chosen by using a modified BIG criterion. 



1 Introduction 



The multivariate statistical theory of shape has been s t udied deeply in the last two decades 

Goodah and Mardia l(|l993[ ). lDi'az-Garci'a eroU (Il997l) iDiaz-Garcia et aZ.I ()2003t ). lDrvden and Mardia 
IQQSf ) and the references there in, Garo-Lopera et all ( 2009() . among many others. Most 



of the works are supported by important restrictions (isotropy) for the covariance ma- 
trix and models (Normal) in order to obtain known polynomial s (zonal polynomials). A 
sort of approache s are g i ven for shape theory, via QR jC oodall and Mardi a (19 9^)), SVD 

LVand KendalTI(ll993l).lGoodalll(ll99ll).lDi'az-Garci^ etalH 200^). afhne (fCoodall and Mardia 

19931 ). iDi'az-Garci'a et all (|2003l ) . ICaro-Lopera et al I (|2009h ). 



Avoiding the restrictions of isotropy and normality carry some problems, because integra- 
ti on over Eucli dean or affine transformations lead to the apparition of invariant polynomials 
of Davis (Il908l) which can not computed for large degrees. 

This work finds a sequence of transformations which let the construction of shape densi- 
ties via the singular value decomposition and based on a non restricted non central and non 
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isotropy elliptical model. The resulting densities avoids the invariant polynomials and they 
are set in terms of series of zonal polynomials which can b e computed by suitable mo difica- 
tions of the existing algorithms for hypergeometric series (jKoev and Edelman 



The work is structured as follows, the main principle and the size and shape distribution 
is given in section[51 then the shape density is obtained in section|31 the associated excluding 
reflection densities are considered in section|3]and finally some particular models are derived 
as corollaries in section [6] which also presents an application in mouse vertebra by studying 
three models (the usual Gaussian and two non Gaussian) with the modified BIG criterion. 



2 Main principle and SVD size-and-shape density 

It is known that the shape of an object is all geometrical information that remains after 
filtering out translation, rotation and scale information of an original figure (represented 
by a matrix X) comprised in N landmarks in K dimensions. So, we say that two figures, 
Xi : N X K and X2 : N x K have the same shape if they are related by a special similarity 
transformation X2 = /3XiH+ljv7', where H : KxK <E SO{K) (the rotation), 7 : Kxl (the 
translation), l^v : N x 1, l^v = (1, 1, ... , 1)', and /3 > (the scale). Thus, in this context, 
the shape of a matrix X is all the geometrical information about X that is invariant under 
Euclidean similarity transformations. 

Now, multivariate statistical theory of shape compares shapes of objects in presence of 
randomness, so if we assume that a figure X, comprised in N landmarks in K, follows 
an elliptical distribution X £7Vx_fs'(Atx' ®! it is of interest to remove translation, 
scaling, rotation from X. Clearly, the sequence LX = Y = H'DP = rW(u)P removes the 
translation (by a sub Helmert matrix L, for example), the rotation (by the SVD of Y) and 
the scale (by dividing for the norm of Y). In order to obtain the density of W we need to 
integrate over the similarity group; it is easy to see that the elliptical assumption lead to 
the product of two tra ces which irre mediably expands in terms of invariant polynomials of 
two matrix arguments ( DavisI ( 19081 )). and the shape densities are not computable for large 



degrees. 

So the classical statistical multivariate analysis restricts the models for the original land- 
mark data in order to obtain densities which are expa nded in terms of studied pol ynomials 
such as the zonal polynomials which are computable (jGoodall and Mardia I (fl993h ). other- 



wise, as we proved in the last sentence, the densities involve non computable polynomials 
for large degrees. 

From the practical point of view the restrictions affect the applications; i.e., the isotropic 
assumption ~ Ik for an elliptical shape model of the form 

X fArxK(/X^, S^, 0, h), 

restricts substantially the correlations of the landmarks in the figure, speciall y in objects with 



symm etries as in the case of mouse vertebra for example, among many others (jPrvden and Mardia 



(|l998l )). So, we expect the non isotropic model, with any positive definite matrix 0, as the 



best model for considering all the possible correlations among the anatomical (geometrical 
o mathematical) points. Howeve r, using the classical approach of the published literature 
of shape (see for example Gooda ll and Mardia ( 1993[ )) under the non isotropic model, we 



obtain immediately invariant polynomials, which can not be computed at this time for large 
degrees. 

In order to avoid this problem consider the following procedure: Let 
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if 0^/^ is the positive definite sq uare root of the matrix 0, i .e. = (0"^^^)^, with 0"^^^ 
K X K, Gupta and Vargal (|l993 . p. 11), and noting that 



where 



then 



X0-1X' = X(0-i/2@-i/2)-ix' = X0-i/2(X0-i/2)' = ZZ', 

Z = X0"i/^ 



with = /x^0"^/^ fsee lGupta and Vargal (|l993l p. 20)). 



And we arrive at the classical starting point in shape theory where the original landmark 
matrix is replaced by Z = X0^^^^. Then we can proceed as usual, removing from Z, 
translation, scale, rotation and/or reflection in order to obtain the shape of Z (or X) via 
the SVD decomposition, for example. 



The SVD decomposition has two version in shape theorv. lGoodalll (|1991I ) and lLe and Kendall 



(|1993[ ). we focus in this paper on Goodall's approach. 

Let n = min(m, K), Y = H'DP be the nonsingular part of the SVD, where H : n x m, 
H e Vn,m represents the Stiefel manifold, D = diag(i?i, D2, . . . , Dn) with Di > D2 > ■ ■ ■ > 
Dn and P : n x K , P e Vn,K- 

Thus the SVD shape coordinates u of X may be found by the following procedure 

LX0^/2 ^ LZ = Y = H'DP = rWP = rW(u)P (1) 

where the SVD shape coordinate system is given by H'D (|Goodalll (|l99lL pp. 296-298)) and 
r = ||HD|| = (trDHH'D)i/2 = (tr 0^)1/2 = ||d|| = ||Y||. Before defining W and u, note 
that when n = K two cases may be distinguished. 

1. P includes reflection, P G 0{k), \P\ = ±1, > and (H,D), written (H,D)^ for 
definiteness, contains reflection SVD shape co-ordinates. 

2. P excludes reflection, P e SO{K), \P\ = +\, \Dk\ > 0, sign{Drn) = sign|X| and 
(H, D), may be written (H, D)^-'^ for definiteness. 

Now the SVD shape matrix W is obtained by dividing the H'D matrix by r, when W may 
include or exclude reflection, in which case we obtain, respectively, = (H'D)-'^/r or 
_ (H'D)^^/r. Finally u is composed of the mn — 1 generalized polar coordinates. 

Our interest now lies in finding the corresponding densities associated with the process 
described in ([T|). Thus we obtain the joint density of (H, D) and the density of W(u). 

In order to obtain the size and shape density we need some integrals involving zonal 
polynomials, extending Jamed ( 1964 . eq. (22)). 



Lemma 2.1. Let X : K x n, Y : K x K and H eVn k- Then 



L 



[tr(Y + XH)]f(HdH') 



Hev„,A- 



/=0 A 



where |(trY)"nrXH| < 1 a^rftrYT^O. 



2. I tr(Y + XH)etr{r(Y + XH)}(HdH') = 

/Hev„,A' 



■etr{rY}{trYoF,C-K::-XX') + J2Y. 

L2"J I ^ * /=0 A ^2 



5) Ca(tXX') 
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where p G 5ft, r € 3?, Ck(B) are the zonal polynomials of ^ corresponding to the partition 
K = (/i, • ■ • /p) of f, with fi = fi a'^^ (a)'^ = Il^=l{a■ ~ U ~ l)/2)/, , (a)/ = a(a + 

1) • • • ( o +.f— 1), are </ie generalized hypergeometric coefficients andgFi is the Bessel function, 
Jame^ ll96A ). 



Proof. 



1. From Lemma 9.5.3 iMuirheadI (jl982l Lemma 9.5.3, p. 397) we have 

/ [tr(Y + XH)]P(HdH') = / [tr(Y + XH)]^'(dH). 

Furthermore, for tr Y 7^ and |(tr Y)-i trXH| < 1 

[tr(Y + XH)]P = (trY)f ^M£(trY)--^(trXH)^. 



Now from Jamed ( 1964 . eqs. (46) and (22))) it follows that 
/ [tr(Y + XH)f (HdH') - 



2-V^»/^g^(,).,(trY)^^/J^^^^^^,^^ 



the result follows, noting that (i)//(2/)! = 1/(4^/!) and that CA(aXX') = a-'^CA(XX'). 

2. This follows b y expanding the exponentials in series of powers and by applying (22) 
and (27) fromjjamed (|l964t) . □ 



Now, the jacobian of the corresponding decomposition is provided next: 

Lemma 2.2. Let be Y : N — 1 x K , then there exist V e Vn,N-i, H e Vn,K and D : nx n, 
D = dia,g{Di,...,Dn), n = min(iV - l),K; Di > D2 > ■■■ > > 6, such that Y = 
V'DH; This factorization is termed non-singular part of the SVD. Then 



(dY) = 2-"|D 



n\Y^\N-l+K-2 



" lliD^ - Dj)idT>)iVdV'){UdU'). 



Proof. See iDiaz-Garcia et ali (|l997f) 



□ 



So, we can obtain: 
Theorem 2.1. The joint density o/(V, D) is 



/v,d(V,D) = 



Tn [f] 

E E ^^TTTTIF^ V'D^V) 



t = K 
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Proof. Let be r2 = S ^fi& ^fi' , so the density of Y is given by 

/y(Y) = -^h [tr (S^^YY' + - 2tr/i'S-iY] . 

|S| 2 

Now, make the change of variables Y = V'DH, so, by Lemma[221 the joint density function 
of V. D. H is 



2-n|j)|JV-l + if-2„ -Q(j)2 _ ^2) 



d^V,D,H(V,D,H) = —j^^ (VdV')(dD) 



xh [tr {T.-^V'T)^V + n) -2tr/i'S"iV'DH] (HdH'). 

Expanding in power series 

dFv.D,H(V,D,H) = (VdV')(dD) 

x^i^*) [tY {-E-^VD^V + n)] [tr(-2tr/x'S-iV'DH)]*(HdH'). 

t=o 



From Lemma \2A\ 

[ [tr (-2trM'S-^V'DH)]^' (HdH') = ^ 7^^^^= (nS-^VD^V) , 

(1) 4* 1 

Observing that = — , the marginal joint density of V, D is given by 

^-lp|^_l+K-2„-Q(^2_^2) 

dFv,D(V,D) = 



rn(f)|s|- 

^ '-^C^n^-'V'-D'V) {VdV'){dB). □ 

t=0 ^-U^/k 

Now, let be 

R' = V'D (2) 

then V = R'D-i, so dV = dR'B^ and VrfV = D-iRdR'D-i. But RdR' : n x n is 

skew symmetric, thus 

(VdV) = |Dr"+i(RdR') (3) 
Theorem 2.2. The SVD reflection size- and- shape density is 

dFRiK) - 2 il^ ^(dD) 

^ ' 7r"(^+^-(3«-i)/2)/2r„ [f ] |S|f 

t=0 ^•(2-^)k 
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Proof. The joint density function of V, D is 
dFv,D(V,D) = 



i<j 



Tn [f\ |S| 



-(dD) 



E E — — ...iJ^ ^ (V'^V) 



Now, let be R' = V'D then V = Il'T) \ so dV = dR'B^ and VdV = D-iRdR'D-i, 
where RdR' is an n x n skew-symmetric matrix. 

(VdV) = |D-i|""^ (RdR') = |Dr"+i (RdR') . 

Thus the joint density function of R, D is 

dFR^D(R,D) = FlFTl^ (d'D) 



r„ [f] 



t=0 



t! (ii^) 



For the integration with respect to D note that if = L, so dLi = 2DidDi and (dL) 
2"|D|(rfD), thus 

(dD) = 2-"|Lr^(dL). 

Let be 



J 



ID 



N+K-3n 



so 



J = j |L3|M-3n-Q(^_^^.)2-«|L|-|(dL) 

= 2-" y |L|^(^+^-3«-i)[](i:,-i:j)(dL) 

From lFang and Zhang I (|l990l ). eq. (3.29), p.l02. we have that 



J = 



- (TV + fc - 2n - i + 1) 



^^n(N+K-2n) + ^n 



(/V + if - 2?i - i + 1) 



+i<'-2n)-n-l]/2 



|L|[(^ 



2-"7r"("^i)/4 JJ^ r 



n(i.-i,)(rfL) 



i=l 



-{N + K -2n-i + l) 



Kjv+g-2«)-S-+ "'"4 



r„ [^(iV + jf-2n)]2-" 
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then 

-p N+K-2n o-n„^ 
^ n 2 



r„ [f]^t(A^+^-^)|s|f 
3 Reflection Shape Density 

For the SVD reflection shape density consider the foUowing transformations 
LX0^/2 = LZ = Y = (V'D)H EE R'H = rWH = rW(u)H, 

where R' = V'D and W = R'/r. 

Now, note that V'D contains {N — l)n coordinates. Then 

vecW = ivec(V'D), r = ||V'D|| = VtiV'D^V = \\Y\\. 
r 



Then, bv lMuirheadI (|l982i Theorem 2.f .3, p. 55) 



(d vec W(u)) = r" [| sin""' d, /\ d9, A dr 

i=l i=l 
m 

= r^'Jiu) /\d9,Adr, 



i=l 



with m = (iV - l)n - f , u = (6*1, . . . , 9m)'- 
Hence, 

Theorem 3.1. The SVD reflection shape density is given by 

dFwCW) = — ^ ^ ' VV ^ , ^ WdW 



X 



/ r"+"+2t-i/j(2t) [^2 WW + tr n] {dr). (5) 

Jo 



Proof. The density of R is 



dF^{R) 



7rt[^+^-'^]r„[f]|s| 



Putting W(u) = R'/r, the joint density of r and u is 

dFrw(u){r,W{u)) = 2_J ^ (r^WdW) 

^ /i(2t) [tr (r^S^i WW + fl)] (r^ffE-'WW) 

t=0 K \ 2 / K 

Note that 
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1. Ck P^s-^ww) = t^'Ck (ns-^w'w). 

2. (r^WdW) = ((rl) WdW (rl)) ^ |rl|"-i (WdW) = r"-i (WdW). 
Collecting powers of r, the marginal of W is 

r„ [ ^+^-^" 1 2-"^^ j(u) , ,^ ^ ^ c« (fis-^ww) 



X 



/ r™+"+2*-i/i(2t) [r2^rS-iWW + trf2] (dr). □ 
Jo 



4 Distributions excluding reflection 

Recall that the SVD shape coordinates u of X are obtained as follows 

LX0^/2 = LZ = Y = HDP = rWP = rW(u)P (6) 



where the SVD shape coordinate system is given by H'D (jGoodalll (|l99ll pp. 296-298)) 
and r = ||HD|| = (trDHH'D)!/^ = (tr 0^)1/2 = ||d|| = ||Y||. When n = K we studied 
the distributions including reflection, i.e. P G 0{k), |P| — ±1, T)k > and (H, D), written 
(H,D)^ for definiteness, contains reflection SVD shape co-ordinates. 

In this section we consider the case when P excludes reflection, thus P G SO{K), 
|P| = +1, \Dk\ > 0, sign(i:)„) = sign|X| and (H,D), may be written (H,D)^-" for 
definiteness. 

Finally, we have that the excluding reflection SVD size-and-shape and SVD shape den- 
sities are given by ^ and ([5]) divided by 2, respectively. 

5 Central Case 

Now, we can derive easily the corresponding central distributions of this work. 
Corollary 5.1. The central reflection SVD size-and-shape density is 

2-np r N+K -2n -\ ^nf 

dFR(R) = — i ^— rr h [trS-^R'Rl (RdR'). 

Proof. Just take /x = in Theorem!^ and recah that h^"^-) = h{-). □ 



Finally, 

Corollary 5.2. The central reflection SVD shape density is invariant under the elliptical 
family and it is given by 

dFw(W) - '"^3 ,f J ^ ^ ftrS-iW'Wl (WdW). 

1/2 

Proof Taking = and s = (trS^^WW) ' r in Theorem 13.11 we obtain the result, 
since 
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^m+n-if^ [r^ trS-^W'W] dr 



= ftrS^^W'W 



h 



m + n 



ds 



(trS^^W'W)- 



6 Some particular models 

Finally, we give explicit shapes densities for some elliptical models. 
The Kotz type I model is given by 



Hy) 





2 






[t 1 + ^^^- 


1)1 



y cxp{-i?y}. 



Then, the corresponding fc-th derivative of h, follows from -j—ry exp{— i??/}, which is 



given by 



dy'^ 



'm — 1 



\{{T-l-^) 



i=0 



see 



Caro-Lopera et 



It is of interest the Gaussian case, i.e. when T — \ and R — j- here the derivation is 
straightforward from the general density. 

/J— s — - 

The required derivative follows easily, it is, h'^^\y) = — kin~i) (~-^)'' 6xp{— and 



r™+"+2t-i/j(2t) p^rS"iW'W + trf2] {dr 



2^^ 

Hence dFw(W) is given by 



■ exp{-i?tr n} (tr S^^WW) 



m + n 



+ t 



^ ,1 .1 N ^ WdW 



t=0 K 



/ ^™+"+2t-i;j(2t) [^2 tr WW + tr n] {dr 
Jo 



i?f-5(™+«)r„ [ ^+^-^" ] J(u) exp{-i?trf2} C„ (flJIS-^WW) 



2n+i7rt(^-^)+^r„ [f ] 
m + n 



x(trs~iww) " r 

Therefore, we have proved that 



+ t 



t=o 

(WdW). 
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CoroUciry 6.1. The Gaussian SVD reflection shape density is 



dFw(W) 



2n+%t(^-^)+f r„ [f ] |S|f 



X (trS-^WW') 
where M = {N - 1)K. 



-si±SL-t „ rm + n 



r[^+t](Wdw'), 



Finally, we propose the result for the Kotz type I model 



h{y) = 







"_fs:(Ar-l)' 

2 




T^K{N- 


-i)/2r 


[t 1 + ^ 


1)1 



y cxp{-i?y}, 



Corollary 6.2. The Kotz type I SVD reflection shape density is 



/w(W) 



-T-H 



[f ] r„ [ ^+^-^" ] J(u)(trJI)^-ietr(-i?f2) 



2«+i7,f (^-^)+f r„ [f ] |s|f r [T - 1 + f ] 

£i: "''""' (w^o(..i:--w-w)- 



t = K 



X < 



u=0 



r[f + f + «]n(r-i-.) 



s=0 



i!i?«(trJ2)»r [T - 1 + f ] 



■E 

m=l 



"m— 1 



n(r-i-i) 



i=0 



(-i?)-"' (tr^2)~"' 
r[T-l-m+f] 



E 

u=0 



r[f + i + J](r-i-m-s) 

s=0 

u!i?"(trfl)« 



where M = {N - 1)K. 

Proof. As we note before the fc-th derivative of h follows from, 

d'' 



dy 



■^y ex.p{-Ry} 



= {-RfyT-^eKp{-Ry}\l+J2 



m=l 



lliT-l-i) 



i-Ry)- 



(7) 



and the corresponding SVD reflection shape density, dFw(W), is obtained after some sim- 
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plification as 



'N+K-2n-\ 



J(u)7r' 



"'0 



2"+%* (^"^)+^r„ [f ] |s|f r [r - 1 + f ] 

.1^1 '-{WdW) (trS-^W'W) 



t = K 



X < 



E 

u=0 



£=0 

!i?"(trO)"r [T - 1 + f-] 



m— 1 



(-i?)-"'(trJ7)" 



E 

M = 



T[T-l-m+f] 
T[f +t + u]Y[{T -1-m- s) 



u-l 



□ 



The Gaussian case can be derived again by taking T = 1 in the above result. 



6.1 Example: Mouse Vertebra 



This classical application is studied in the Gaussian case bv iDrvden and Mardia I (|1998[ ). 
Here we consider again the same model and contrasted it, via the modified BIG criterion, 
with two non Gaussian models. 

The isotropic Gaussian shape density is given by 



dFw W = — n 2 i ^ / ctr - — 

^ ' ^^{2M+n-3n^+2nM)^-(m-M+n)Y^^K^ \ 2a 



2 



X 



t = K \2^>K 



(WdW), (8) 



where M — K{N — 1), n = min{(A^ — 1), K} and m — [N — l)n — 1. Here we study three 
models, the Gaussian shape (N), and the Kotz (K) model for T = 2 and T — i. 

The shape density associated to the Kotz model indexed by T = 2, R — i (and s = 1) 
is given by: 
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dFw{W) 



■ etr 



(tr (f^) - 2t) r + <] + r [2i±2i + 1 + 1] 



t=0 

K V 2 / At 

And the corresponding density, dFw(W), for the Kotz model T = 3, is obtained as: 



^i(2M+«-3«2+2„M)^-(™-M+n) 2)r„ [f " 



■ etr - 



2(72 



E 



4r - 2t - 4itr 



-At + 2 tr 



2(72 



V2(72 

m + n 



V2(72 

i + 1 



TO + n 



2 

TO + n 



+ t + 2 



^ ^2^^^^^^ tfZ(W(iW'). 



(9) 



In order to decide which the ehiptical model is the best one, different criteria have been 
employed fo r the model selection. We shall consider a modification o f the BIC statistic as 
discussed in lYang and Yand ( 2007 ) , and which was first achieved by iRissane m a 

coding theory framework. The modified BIC is given by: 

BIC* = -2£(/I, ct2, /i) + np(log(n + 2) - log 24), 

where ?2, K) is the maximum of the log- likelihood function, n is the sample size and Up 
is the number of pa rameters to be estimated for ea ch particular sh ape density. 

As proposed by iKass and Rafterv ( 19951 ) and Rafterv ( 1995 ). the following selection 
criteria have been employed for the model selection. 



Table 1: Grades of evidence corresponding to values of the BIC* difference. 



BIC* difference 


Evidence 


0-2 


Weak 


2-6 


Positive 


6-10 


Strong 


> 10 


Very strong 



Fixing the variance of the process as 50 (the maximum median variances of the two 
samples), the maximum likelihood estimators for location parameters associated with the 
small and large groups are summarized in the following table: 

According to the modified BIC criterion, the Kotz model with parameters T — i, R — ^ 
and s = 1 is the most appropriate among the three elliptical densities for modeling the data. 
There is a very strong difference between the non Gaussian and the classical Gaussian model 
in this experiment. 
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Table 2: The maximum likelihood estimators 



Group 


BIC* 


Mil 


Ml2 


M21 


At22 


M31 




a 

K ■T = 3 














if :T=2 














-5.9146 


-1.9214 


-42.5338 


14.1761 


-4.8190 


24.0766 


Small 


-39.6272 


2.1250 


-47.8016 


16.3513 


-3.9691 


26.9744 




-23.0250 


-6.1682 


-45.0331 


14.6035 


-6.4983 


25.5710 




24.5000 


-23.1834 


-32.9246 


8.1455 


-10.5612 


22.5117 


Large 


-9.2156 


-34.4190 


-29.3792 


5.9227 


-13.7687 


20.9197 




7.3880 


-42.9562 


-1.6455 


-3.6885 


-13.7522 


3.6578 



M32 


^41 


M42 


M51 


M52 


-0.2605 

2.1288 

-2.6121 


4.8195 

4.9123 

5.5892 


4.7981 

5.8552 

4.6337 


-29.2691 

-33.0646 

-30.8241 


-0.7072 

0.4382 

6.0498 


-12.8674 
-20.2846 
-27.4382 


5.0300 

6.1263 

5.2368 


2.3324 

1.1366 

-2.7585 


-26.0251 
-23.1221 

-0.9942 


18.6865 

27.5996 

34.1886 



Let fii and be the mean shape of the small and large groups, respectively. We test 
equal mean shape under the best model, and the likelihood ratio (based on —2 log A « Xio) 
for the test Hq . fii = ■ Mi 1^2^ provides the p- value 0.84, which means that there 

are extremely evidence that the mean shapes of the two groups are equal if the variance of 
the experiment is fixed in 50 (the maximum median of the variances of the two samples), a 
deeper study of this case is suggested, because the variance estimation was problematic in the 
performed inference procedure for these data. We highlight that our intention is to illustrate 
the technique and performed inference with an exact likelihoo d efficiently computable aft er 



modification of the algorithms given for hypergeometric series (jKoev and Edelman l (|2006l) ). 

A final comment, for any elliptical model we can obtain the SVD reflection model, 
however a nontrivial problem appears, the 2i-th derivative of the generator model, which 
can be seen as a partition theory problem. For The general case of a Kotz model (s ^ 1), and 
another models like Pearson II and VII, Bessel, Jensen-logistic, we can use formulae for these 



derivatives given by ICaro-Lopera et al ( 2009t ) . The resulting densities have again a form 



of a generalized series of zonal polynomials which can be co mputed efficiently af t er som e 



modification of existing works for hypergeometric series (see Koev and Edelman ( 2006l )'l. 



thus the inference over an exact density can be performed, avoiding the use of a ny asymptoti c 
distribution, and the initial transformation avoids the invariant polynomials of lDavisl(ll908l) . 
and it lets the inclusion of any correlation among landmarks. 
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